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I. INTRODUCTION 

The KPZ problem [l| describing the growth in time of an interface in the presence of noise have been the subject 
of intense investigations during the past almost three decades dSl. During last five years the major progress in 
solving this problem has been achieved which resulted in derivation of the exact solutions for various types of the free 
energy probability distribution function (PDF) 

In this paper I would like to consider the possibility to apply the ideas and technique developed in the studies 
of the KPZ problem for formally equivalent problem of the randomly forced Burgers turbulence (the ’’Burgulence” 
problem). Here one considers a velocity field v{x,t) governed by the Burgers equation 

dtv(x, t) + v{x, t)dxv(x, t) = vdlv{x, t) + ({x, t) (1) 

where the parameter i/ is the viscosity and C,{x,t) is the Gaussian distributed random force which is J-correlated 
in time and which is characterized by finite correlation length ^ in space: C,{x,t)C,{x',t') = u5{t — t')U\{x — x')/^]. 
Here U(x) is a smooth function decaying to zero fast enough at large arguments and the parameter u is the injected 
energy density. In this problem one would like to derive e.g. the the probability distribution functions of the velocity 
gradients P[dxv{x,t)] or two-points distribution function P[v{x,t),v{x',t)] at scales smaller than the length scale ^ 
of the stirring force C (see e.g. and references there in). 

Formally the above problem is equivalent to the KPZ equation as well as to the (1+1) directed polymers. Indeed, 
redehning v{x,t) = —dxF{x,t) and Q{x,t) = —dxV{x,t) one gets the KPZ equation for the interface profile F{x,t) 
(which is the free energy of (1+1) directed polymers): 

dtF{x,t) = ^{dxF{x,t)y - ^Td^F(x,t) - V(x,t) (2) 

where T = 2^ is the temperature parameter of the directed polymer problem and V(x, t) is the Gaussian distributed 
random potential. 

The idea of a new approach to the Burgulence problem which I would like to demonstrate in this paper is in the 
following. According to the above definitions the velocity field can be represented as 

v{x, t) = hm[F'(x + e, t) — F{x, t)]/e (3) 

e—>-0 

Thus, deriving the four-point KPZ probability distribution function V[F{x + e,t), F{x,t), F{x' + e,t'), F{x',t')] and 
taking the limit e —>• 0 one could hopefully obtain the result for P[v{x,t),v{x' The only ’’little problem” is that 
unlike the usual KPZ studies operating with the (^-correlated in space random potential, in the Burgulence problem one 
is mainly interested in the spatial scales comparable or much smaller than the random potential correlation length 
In other words, in this approach, first one has to study KPZ problem with random potentials having finite correlation 
length. In the present study (as a matter of ’’warming up” exercise) I’m going to consider another ’’extreme case” 
in which the random potential V{x,t) of the KPZ problem is changed by it’s linear approximation: V{x,t) —> C (t)x 
where ^(t) is Gaussian distributed random force. In this case we obtain the model introduced by Larkin |37l . l38l | 
long time ago to study small scale displacements of directed polymers. In this approximation the model becomes 
Gaussian and therefore exactly solvable. Nevertheless, the statistical properties of its free energy (as well as some 
others quantities) turn out to be rather non-trivial (see e.g (39l - l^ ). For that reason this model hopefully could serve 
as a good ground for testing various approaches developed in the recent KPZ studies. 
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In Section II we introduce the model and formulate the main ideas of the replica approach which will be used in 
the further calculations of the two-time free energy distribution functions (for the directed polymer model) and the 
corresponding velocity distribution function (in the Burgulence problem). In Section III as the matter of the demon¬ 
stration of the replica technique the explicit expression for the two-time (two-point) free energy distribution function 
is derived. In Section IV the two-time four-point free energy distribution function is calculated and the corresponding 
two-time velocity distribution function P[v{x, t), v{x',t')] of the Burgulence problem is obtained. Further perspectives 
of the present approach is discussed in Section V. 


II. THE MODEL AND THE METHOD 


In this paper we consider the model of one-dimensional directed polymers defined in terms of an elastic string 4’{t) 
directed along the r-axes within an interval [0, t] which passes through a random medium described by a random 
potential V{4i, r) = ((T)(j). The energy of a given polymer’s trajectory (/>(t) is given by the Hamiltonian 


H[(j){T),V] = f dT^^^[^r(|){T)]'^ + CiTM'r)]; 


( 4 ) 


where the random force C{t) is described by the Gaussian distribution with a zero mean C(t) =0 and the d-correlations: 

= uS{t - t') (5) 

The parameter u describes the strength of the disorder. 

For the fixed boundary conditions, (()(0) = y, = x, the partition function of the model Q is 


r4>{'t)=x 

Z{x\y\t)= / D(/)(t) = exp[-l3F{x\y,t)] 

d<l>{0)=y 


( 6 ) 


where /3 is the inverse temperature and F{x\y;t) is the free energy. In the replica approach one consider the average 
of the iV-th power of the above partition function: 


Z^{x\y,t) = Z{N;x\y]t) = exp[-l3N F{x\y]t)] 
where (...) denotes the average over the random force Simple Gaussian averaging yields: 


N 


Z{N-,x\y,t) = 


a=l 


r4'a{t)=x 



/ T>())a(T) 

exp 


d^a{0)=y 




where 


1 /■* / ^ AT 

Hn[(I}] = 7: rfT y'[^r</>a('r)]^ - /3u V((>a(T)(;i6(r) 

•’'o \a=l a,b 


( 7 ) 


( 8 ) 


( 9 ) 


is the Gaussian replica Hamiltonian. 

Introducing the free energy distribution function, Px\y,t{F), the relation ([7]) can be represented as follows: 

/ + 00 

dP Px\y,t{F)exp[-l3NF] (10) 

-OO 

which is the Laplace transform of the distribution function, Px\y,t{F) with respect to the parameter (dN. In the 
lucky case when the moments of the partition function Z[N]x\y,t) allows an analytic continuation from integer to 
arbitrary complex values of the replica parameter N the above relation makes possible to reconstruct the probability 
distribution function Px\y-t{F) via the inverse Laplace transform: 

/■+*“ s 

Px\y,t[F) = J —z{--x\y,t) exp(sF) (11) 

In the present model this distribution can be computed explicitly and the resulting function Px\y,t{F) turns out to be 
rather non-trivial (sol - l^ . 
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Besides the replica partition function eq.® it is convenient to introduce the A^-particle ’’wave function”: 

N 


«'Ar[x|y;t] = 





/ V4>a{x) 

y<l>a.{0) = ya. 

exp 

-fiHN[(p'] 


( 12 ) 


where the Hamiltonian Hn[4>\ is given in eq.®. In the present model the function 'I'Ar[x |y;^] can be computed 
explicitly (see e.g. (4ll|l: 


N 


^Ar[x|y;t] = C{N,t) exp<j - VaY - ^A{N,t) 


where 


2t 

C{N,t) 

A{N,t) 

B{N,t) 


a—1 


2t 


N 


J2{Xa-ya) 


.a—1 


N 


+ jB{N,t) ^ XaVb} (13) 


,,b=l 




\2TTtJ sin(y^ /SNuf^) 

1 / _ \ 
N y tan(-\/ PNut^) j 


1 [j _ cos(^/:8^^)] 

sin(Y^ l3Nut‘^) 


(14) 

(15) 

(16) 


Note that the above expression for the wave function ea. dT^ is valid only at finite time interval: t < = 

tc{N). The reason is that dne to specific form of the interaction potentials in the replica Hamiltonian ([9]) the 
directed polymers trajectories go to infinity at finite time t r(N). For the original physical system this phenomenon 
is explained by the presence of the slowly decaying left tail [39l - l4ll| of the free energy distribution function P^\y.t{F) 
which (according to the relation (I10|) l results in the divergence of all partition function moments Z{N; x\y;t) with 
N > N,{t) = ^(/3nt2)-i 


III. TWO-TIME FREE ENERGY DISTRIBUTION FUNCTION 


For simplicity, in this section we consider the directed polymer problem with the zero boundary conditions: <('(0) = 
(j){t) — 0. For a given realization of the random function (((r) let us denote by Fi the free energy of the directed 
polymers with the (zero) ending point at time t and by F 2 the one of the directed polymers with the (zero) ending 
point at time t + At. According to the definition, eq.dH]), for the difference of these free energies, F = F 2 — Fi one 
has: 


exp{-/3F’} = Z-i(0|0; t)Z(0|0; t +At) 


(17) 


Taking n-th power of the above relation and performing the averaging over quenched randomness we gets 

exp{-/3nF’} = Z-"(0|0; t) Z"(0|0; t + At) (18) 

Introducing the two-time free energy difference probability distribution function in terms of the replica 

approach the above relation can be represented as follows: 


^ + 00 


dF Pt^At{F) exp{-finF} 


hrn^ Z^-"(0|0; t) Z^(0\0; t + At) 


(19) 


where according to the usual replica formalism, in the first step, the r.h.s of the above relation is computed for an 
arbitrary integer N and then, at the second step, the obtained resnlt is analytically continued for arbitrary real N, 
and finally the limit TV —)> 0 is taken. 

For an integer N > n m terms of the wave function, eqs. (IlU-dlSl), the prodnct of the partition functions in the 
r.h.s. of ea. dTOl) can be represented as follows: 


Z^-"(0|0; t)Z"(0|0; t + At) = Z{N, n-1, At) 


/ +00 

dxi...dXn'i>N\^--,0,Xi, ...,Xr,\0;t] 4'* [xi, ...,x„|0; At] 

-00 


( 20 ) 


N—n 
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FIG. 1: Schematic representation of the directed polymer paths corresponding to eq. (O 


where the second (conjugate) wave function represent the ’’backward” propagation from the time moment {t + At) to 
the previous time moment t. Schematically the above expression is represented in Figure 1. 

Substituting ea. lfT^ into ea. lQlH) we get: 

/ + 00 ( 1 pi I Ai 1 1 1 'I 

dxi...dx„exp<--P ^ +-A(iV,t) + —yl(n, At) XaXt> (21) 

I 0.6=1^ J 

Simple integration yields: 


Z{N, n-1, At) = C{N, t)Cin, At) 


2-KtAt l"/2 


1 + n 


tAt /I 


1 




L I3{t + Ai) - 

Substituting here the explicit expressions (m and (USD, and taking the limit A — 0 we obtain 

rt/2 


- 1/2 


lim Z{N ,n]t, At) = ( -—) 

N^O ^ ^ \t + AtJ 

substituting ea. (l23l) into ea. dT^ and redefining: 


s/l3nu{At)'^ (t + At) 


\ sin(-y/3nii(At)2) At — \l3nu{At)t‘^ + t — 

\ VV \ ; L 3^' V 1 tan(^/3nu(At)2) 


At = (t 
Pnu{At)^ = u! 

F = uet^f 

we get the following relation for the probability distribution function Pt,^{f) of the rescaled free energy /: 


/ +00 ^ 

d/Pt,a/)exp{-cc/} = (1 + e)"^ 

-OO 




sin(v^) + ^/u} cos{y/Lo) 


( 22 ) 

(23) 

(24) 

(25) 

(26) 

(27) 


By inverse Laplace transform in the limit when both t —^ oo and At —>■ oo (such that the parameter (^ = At/t remains 
finite) we get the following universal result for the limiting two-time free energy distribution function: 


/ +io 

-zoo 


27rz 


a;i/4exp{w/} 


\l sin(Vw) + -v/^cos(z/w) 


( 28 ) 
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It is interesting to note that this function (like its one-time counterpart [39l - l4l| ') is identically equal to zero at / > 0. 
Indeed, since at / > 0 the function under the integral in the r.h.s of ea. (E51) quickly goes to zero at w —>■ —oo, the 
contour of integration in the complex plane can be safely shifted to —oo, which means that > 0) = 0. 


IV. TWO-TIME VELOCITY DISTRIBUTION FUNCTION 


Velocity in the Burgers problem is given by the derivative of the free energy of the directed polymer problem: 

, dF{x,t) F{x + e)- F{x,t) 

v(x,t) = --- = — hm- (29) 

ox e->-o e 

Thus, to compute the two-point velocity distribution function in terms of the directed polymers, first, keeping e finite 
we consider specially constructed four-point object (see below), and only in the final stage of calculations we take the 
limit e —>■ 0. 

According to the relation ([6]) 


exTp{—l3[F{x + e,t) — F{x,t)]} = exp{l3ev{x,t)} = 


Following the procedure described in the previous section we have: 


Z{x + e|0; t) 
Z{x\()-, t) 


exp{/3niet!(a;i, t) -|- I3n2ev{x2,t + At)} = 


(30) 


(31) 


= lim Z"i(xi + elO; t) Z^i“"i(a:i|0; t)Z^‘^(x 2 + e|0; t + At) Z^'^~'^^(x2\0- t + At) 

JVi,Af2->0 

Introducing two-time velocity distribution function Pxi,x 2 t,Ativi,V 2 ) the above relation can be represented as follows: 

f*+oo 


where 


dvidv 2 Pxi,x 2 t,Ati'>^i^'^" 2 )exp{Pnievi+ I3n2ev2} = lim Ze(Vi, m, V 2 ,ni,xi, 0 : 2 , t, At) (32) 

^1 ,A^2~^0 


Z^[Ni,ni,N2,n2,xi,X2,t,At) = (33) 


= Z”i(a;i -I- e|0; t) Z^^~'^'^(xi\0\ t)Z^^{x 2 + e|0; t -|- At) Z^‘^~'^'^{x2\0\ t + At) 
In terms of the wave function, eq. ©-([ISl), 

Z<:(Ai,ni, V2,n2,a:i,a:2,t, At) = 

/ +00 r+co 

Dy / Dz '^Ni+N 2 ■■■,xi, xi + e, . . .,xi -I- e ,yi, ■.■,yN 2 — 712 1 ^712 I ^ 

-OO ^—00 ^ V ^ V 

Ni—Tli Til 


where 


'X'^*N2[yi’ ■■■’yN2-n2, Zl, Zn^ \ X 2 , . . ^.,X2 , X2 + 6, ..^.,X2 + C ; At] 

N 2 — 7 I 2 '^2 


f + OO N2—7I2 p~\-oo 

Dy = n / 

-00 d -00 


^ + 00 


D. s n 


^2 /*+oo 


dZa 


(34) 


( 35 ) 


Schematically the expression in ea. (l34l) is represented in Figure 2. 
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FIG. 2: Schematic representation of the directed polymer paths corresponding to eq. dSl 


Substituting the explicit expressions for the wave function m into ea. (l34l) we get: 


Ze{Ni,ni,N2,n2,Xi,X2,t,At) 


C{Ni + N 2 , t) C{N2, At) X 


r+oQ p+oo I p 

X / Dy / Dz exp< — — 


' —00 


N2—n2 n2 


-AA{N,+N2,t) 


P 'N2—n2 

'2Ai 


{Ni-ni)xl+ni{xi+ef + ^ yl + ^z, 

a=l a—1 

N2 — n2 n2 

{Ni — ni)xi + ni{xi + e) + E 
n2 

iVa - X2f + -X 2 -e) 

a—1 

12—^2 n2 

Y “ ^ 2 ) + E(^“ -X2-e) 

a=l a—1 

„ N2 — n2 n2 'j 

+ ^S(A2,At) (E 2^“ + E^“)[ {N 2 — n 2 )x 2 + n 2 {x 2 + e)j > 

a—1 a—1 ) 


a=l 


(36) 


where A{N,t), B{N,t) and C{N,t) are given in egs. (ITU) - ([TO . Introducing fV 2 -component vector x = 
{yi, ■■■,yN2-n2,zi, Zn^} after simple algebra we get 

Z£(fVi,ni,iV 2 ,n 2 ,a;i,a; 2 ,t, At) = C{Ni + N 2 ,t) C{N 2 , At) x 


where 



^ N 2 N 2 

2 'y ^ TabXaXb d" ^ ) L^Xa 

a,h—l a=l 



(37) 


G = 


+ 


i [ni(xi + e)^ + (tVi - ni)xl] 
jA{Ni + N2,t){NiXi + nie)^ 


^ [n2(x2 + e)^ + (A 2 - n2)x2] + 
■^A{N2,At){N2X2 + n2e)^ 


( 38 ) 
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and 


Tab = 


1 = P 


P 


iSab + K 

t + At 


tAt 


p 


K = ^A{N,+N2,t) + ^A(7V2,A0 


La — L + Xa 


(39) 

(40) 

(41) 

(42) 


L = —^A{Ni + N2,t){LfiXi + nie) + — [A(iV2, A<) + 5(^2, Al)] (iV2X2 + n2e) 


= 


X 2 + e, for 0=1,..., 77,2 


X2, for O = 772 + 1, ■■■,N 2 
Simple integration over x’s in eq. yields: 

Z^[Ni,ni,N 2 ,n 2 ,xi,X 2 ,t,^t) = C{Ni + N 2 ,t) C{N 2 , At) x 


Nz 


1 1 - 1 \ 
X exp-^--/JG--Tr InT + ^ X! 


ab 


,,b^l 


where 




Tr In T = N 2 In 7 + In (1 + N 2 


Tab" = -5ab - 


7 ° + N2k) 

Simple calculations yield: 

Z^[Ni,ni,N 2 ,n 2 ,xi,X 2 ,t,^t) = G(Ai + A2, t) G(A2, At) x 


X exp|-i/3G- ilV2ln7- iln(^l + ^2^) 


L^N- 


2 , PL 


2(7 + N 2 H) jAt 

PkN2{N2X2 + 7726 ) P‘^k{N2X2 + 7726 )^ / 3 ^ 


+ 71 X 7 (^ 20:2 +7726)- 


{N 2 - 772 ) 0:2 + 772 ( 0:2 + e)" 


7At(7 + A 2 K) 27(At)2(7 + A^2'«) 2j{At)‘^ . 

Next step of the calculations is to take the limits N 12 —>■ 0. Using explicit expressions PX1) -P ^ . 
, one easily finds: 

lim C(N, t) = 1 

7V->0 

lim A(N, t) = - But^ 

N^O 3 

lim B{N,t) = -But^ 

TV^o ^ ’ 2 

lim [BG = ^ " (2xi+e) + '^.^ (20:2 + e)- {Pnie)‘^ut - {pn 2 e)‘^uAt 

Ni,N2^o t At 3 3 

lim L = -B^meut + -P^n 2 euAt 

Ni,N2^o 3 6 

lim K = —-/3^u(t + At) 

Ni,N2^o 3 


(43) 

(44) 


(45) 

(46) 

(47) 


(48) 
, dH]) and 

(49) 

(50) 

(51) 

(52) 

(53) 


( 54 ) 













Substituting the above limiting values into ea. (|48)) we get 


lim Z^[Ni,ni,N 2 ,n 2 ,Xi,X 2 ^t,M) = Zf^(ni,n 2 ,Xi,X 2 ,t, At) = 

A^1,2^0 

= + liPn.erut + liPn 2 e)Mt + At) + ( 55 ) 

Substituting the above result into ea. (l5^ and introducing notations /3ni_2e = Si ,2 in the limit e —>■ 0 we obtain: 


+ 00 


dvidv 2 Pxi,x 2 t,Ativi,V 2 ) exp{sivi + 52 ^ 2 } = 


Xl 

= exp-; -—si - 


X 2 1 n 1 , , , 9 Ut^ 

- 7-S2 + -xutsi + -uit + At)s2 + 777-T-^ 

t + At 6 ^ 6 ^ ^ 3(t + At) 


S 1 S 2 


Redefining 


Si = \ — CVi 

ut 


S2 = 


u(t + At) 


Ui = — 


Xl 


UJ2 


ut Vi 


V2 = - 


X2 


t At 


+ \ -u(t + At) V2 


At = ^t 


(56) 

(57) 

(58) 

(59) 

(60) 
(61) 


we get the following relation for the probability distribution function 'P^{vi, 'O 2 ) for the rescaled velocities vi and V 2 , 

eas.(l5^-(l5ni): 


f /■+°° r 1 1 

J J dviV2'P^{vi, V2) exp{ujivi + u;2V2} =exp|-wj + - 


2 , wia;2 
LO 2 


Performing simple inverse Laplace transformation 


V^iyi, 'D 2 ) = 


+ ZOO 


diuidiU2 


>{) 


1 2 I f 2 , W 1 W 2 


(l+^)3/2 


(l+^)3/2. 


■ UJiVi - UJ2V2 


c (27rt)2 

one eventually obtain the following very simple result for the two-time velocities distribution function: 

1 / (1+0' „ f (1+0' 


V 2 ) = 7 ^ 


2 ^ Y (1 + ^)3 _ 1 2 [(1 + - 1 ] ^ (1 + 0 '/" 


(«? 


- 2 - 


V 1 V 2 


+2 


(62) 


(63) 


(64) 


where ^ = At/t is the reduced separation time parameter. 

One can easily check that in the limit of infinite separation time, ^ > 00 , the distributions of two velocities are 

getting independent: 


YimP,(fii,fi2) = ^ exp{-ifi?-ifi2^} 

while in the opposite limit of coinciding times, ^ ^ 0, one finds 

^2) = exp|-ifii| i 5 (-Di -V2) 


(65) 


( 66 ) 
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as it should be. 

Besides, using the exact result, ea. dMl) one can easily compute the time dependence of the two velocities correlation 
function: 


as well as the probability distribution function for the velocities difference v = V 2 — vi: 

2^ Y47r[(l + 03/2-1] '"’'I 4[(1 + 03/2-1] ' 

V. CONCLUSIONS 

In this paper we have considered the problem of velocity distribution functions in the Burgulence problem in terms 
of the toy Gaussian model of (1+1) directed polymers. In particular the exact result for the two-time free energy, 
eg. ( 1 ^ 51 ) . and two-time velocity distribution functions, ea. dMl) has been derived. Of course the considered system is too 
far from the realistic one. Nevertheless, it has one important advantage: being exactly solvable, some of its statistical 
properties are rather non-trivial. All that, in my view, makes this model to be rather useful tool for testing new ideas 
and various technical aspects of the calculations (like the replica technique considered in this paper). 

Following the proposed route, the next step would be to consider the model with finite range correlations of the 
random potentials. Of course, one can not hope to get exact results here. In terms of the replica approach, first of 
all, one is facing the problem of A^-particle quantum bosons with attractive finite range interactions whose solution 
is not known. Nevertheless even the qualitative understanding of the structure of the N-particle wave function of 
this system (which at the qualitative level might be not so much different from that of the Bethe ansatz solution for 
the (5-correlated potentials) could hopefully be sufficient to get some understanding of the velocity statistics in the 
Burgulence problem. 


(67) 

( 68 ) 
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